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Abstract 

Parity-violating (PV) interactions among quarks in the nucleon induce 
a PV ^yNN coupling, or anapole moment (AM). We compute electroweak 
gauge- independent contributions to the AM through 0(1/A^) in chiral per- 
turbation theory. We estimate short-distance PV effects using resonance 
saturation. The AM contributions to PV electron-proton scattering slightly 
enhance the axial vector radiative corrections, R P A , over the scale implied by 
the Standard Model when weak quark-quark interactions are neglected. We 
estimate the theoretical uncertainty associated with the AM contributions 
to R A to be large, and discuss the implications for the interpretation PV of 
ep scattering. 

PACS Indices: 21.30.+y, 13.4fl.Ks, 13.88.+e, 11.30.Er 



I. INTRODUCTION 



The SAMPLE collaboration at MIT-Bates has recently reported a value for the 
strange-quark magnetic form factor measured using backward angle parity-violating (PV) 
electron-proton scattering 

G^(Q 2 = 0.1 GeV 2 /c 2 ) = 0.61 ±0.27 ±0.19 , (1) 

where the first error is experimental and the second is theoretical. The dominant contri- 
bution to the theoretical error is uncertainty associated with radiative corrections to the 
axial vector term in the backward angle left-right asymmetry A LR [[J: 

A LR cc Q p w + QZ^- + Q|?^S- - (1 - 4sin 2 9 W ) y/l + , (2) 

where and are the proton and neutron weak charges, respectively, Q$ is the 
SU(3)-singlet weak chargeQ, 9 W is the weak mixing angle, and r = Q 2 /4M 2 . The axial 
form factor is normalized at the photon point as 

G>M = -9A[i + nZ\ (3) 

where g& = 1.267 ±0.004 || is the nucleon's axial charge as measured in neutron /3-decay 
and EP A denotes process-dependent electroweak radiative corrections to the V(e) x A(p) 
scattering amplitude. 

The radiative correction R A is the subject of the present study. It was first analyzed 
in Ref. and found to be large, negative in sign, and plagued by considerable theoretical 
uncertainty. Generally, R P A contains two classes of contributions. The first involve elec- 
troweak radiative corrections to the elementary V(e) xA(q) amplitudes, where q is any one 
of the quarks in the proton. These radiative corrections, referred to henceforth as "one- 
quark" radiative corrections, are calculable in the Standard Model. They contain little 
theoretical uncertainty apart from the gentle variation with Higgs mass and long-distance 
QCD effects involving light-quark loops in the Z — 7 mixing tensor. The one-quark con- 
tributions can be large, due to the absence from loops of the small (1 — 4sin 2 6 W ) factor 
appearing at tree level (see Eq. (^) ) and the presence of large logarithms of the type 
\n(m q /M z ). 

A second class of radiative corrections, which we refer to as "many-quark" corrections, 
involve weak interactions among quarks in the proton. In this paper, we focus on those 
many-quark corrections which generate an axial vector coupling of the photon to the 
proton (see Figure 1). This axial vector ppy interaction, also known as the anapole 
moment (AM), has the form 



*Note that in Ref. (2|, the weak charges are denoted £ 
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c am = ±-N(a s + a v r 3 ) W Nd u F^ . (4) 

x 

(Here, we have elected to normalize the interaction to the scale of chiral symmetry break- 
ing, A x = 4TrF n .) These many-quark anapole contributions to R A , which are independent 
of the electroweak gauge parameter 0, were first studied in Ref. ||,|5| and found in Ref. 
H] to carry significant theoretical uncertainty. The scale of this uncertainty was estimated 
in Ref. |J, and this value was used to obtain the theoretical error in Eq. (H). (Note that 
the central value for Gm given in Eq. ([!]) is obtained from the experimental asymmetry 
using the calculation of Ref. ||). 

In order to better constrain the error in G$ associated with R A1 the SAMPLE col- 
laboration performed a second backward angle PV measurement using quasielastic (QE) 
scattering from the deuteron. The asymmetry A ifl (QE) is significantly less sensitive to 
Gm than is A LR (ep), but retains a strong dependence on R^ =1 , the isovector part of R P A . 
The calculation of Ref. || found the uncertainty in R A to be dominated by this isovector 
component — R^ =1 ~ —0.34 ± 0.20 — and the goal of the deuterium measurement was, 
therefore, to constrain the size of this largest term. A preliminary deuterium result was 
reported at the recent Bates25 Symposium at MIT, and suggests that R A =1 has the same 
negative sign as computed in Ref. [Q but has considerably larger magnitude, possibly 
of order unity ||. Combining this result with the previous A LR (ep) measurement would 

(s) 

yield a nearly vanishing value for G M , rather than the large and positive value quoted in 
Eq. f). 

The prospective SAMPLE result for -R^ -1 is remarkable, indicating that a higher-order 
electroweak radiative correction is of the same magnitude as, and cancels against, the tree- 
level amplitude! The occurance of such enhanced electroweak radiative corrections is rare. 
Nevertheless, there does exist at least one other instance in which higher-order electroweak 
processes can dominate the axial vector hadronic response, namely, the nuclear anapole 
moment. The anapole moment of a heavy nucleus grows as 

A 2/3 

(see, e.g. Refs. |§,[j],[§ and 

references therein). Because of the scaling with mass number, the nuclear AM contribution 
to a V(e) x v4(nucleus) amplitude can be considerably larger than the corresponding tree- 
level Z°-exchange amplitude, and this A 2 ^ 3 enhancement is consistent with the size of 
the Cesium AM recently determined by the Boulder group using atomic parity-violation 
|10|| . The reason behind the enhancement of R^ =1 for the few-nucleon system, however, 
is not understood. The goal of the present paper is to investigate whether there exist 
c onventional, hadronic physics effects which can explain the enhancement apparently 
implied by the SAMPLE deuterium measurement. 



In order to address this question, we revisit the analysis of Ref. 0. Following Ref. [TT], 
we re-cast that analysis into the framework of heavy baryon chiral perturbation theory 
(HBChPT) 013. We carry out a complete calculation of R T A =l and R T A =0 to order 1/AJ, 



including loop diagrams not considered in Refs. . We also extend those analyses to 
include decuplet as well as octet intermediate states, magnetic insertions, and SU(3) chiral 
symmetry. As in Ref. ||, we estimate the chiral counterterms at (9(1/A^) using vector 
meson saturation. However, we go beyond that previous analysis and determine the sign 
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of this vector meson contribution phenomenologically. We find that decuplet intermediate 
states and magnetic insertions do not contribute up to the chiral order at which we trun- 
cate. Also, the effect of SU(3) symmetry, in the guise of kaon loops, is generally smaller 
than the pion loops considered previously. In the end, we express our results in terms of 
effective PV hadronic couplings. Some of these couplings may be determined from nu- 
clear and hadronic PV experiments or detailed calculations (for reviews, see Refs. flI3UI5| ), 
while others are presently unconstrained by measurement. Guided by phenomenology and 
the dimensional analysis of Ref . |TTfl , we estimate the range of possible values for the new 



couplings. We suspect that our estimates are overly generous. Nevertheless, we find that - 
even under liberal assumptions - the AM contributions to R^ =l appear unable to enhance 
the one-quark corrections to the level apparently observed by the SAMPLE collaboration 
and, in our conclusions, we speculate on possible additional sources of enhancement not 
considered here. 

The remainder of the paper is organized as follows. In Section 2, we relate the 
anapole couplings a S)V to the radiative corrections, R^ =0,1 , and in Section 3, we outline our 
formalism for computing these couplings in HBChPT. A reader already familiar with this 
formalism may wish to skip to Section 4, where we compute the chiral loop contributions to 
the nucleon anapole moment through 0(1/A?.). We also include the leading l/m N terms 
in the heavy baryon expansion, which generate contributions of 0(l/A x m N ). Section 
5 contains the vector meson estimate of the chiral counterterms and the determination 
of the sign, while Section 6 gives our numerical estimate of the AM contributions to 
rt=o,x ^y e k r i e fly discuss the phenomenology of hadronic and nuclear PV and what that 
phenomenology may imply about the scale of the unknown low-energy constants. Section 
7 summarizes our conclusions. The Appendices give a detailed discussion of (A) our 
formalism, (B) the full set of hadronic PV Lagrangians allowed under SU(3) symmetry, 
and (C) graphs, nominally present at 0(1/ A£) but whose contributions vanish. 



II. ANAPOLE CONTRIBUTIONS TO R A 

The electron-nucleon parity violating amplitude is generated by the diagrams in Figure 
2. At tree level this amplitude reads 



where 



and 



iM pv = iMZ + iMZ , (5) 



iM% = i^j=l>*< N\J X \N> (6) 



tMZ = t^l X < N\J X5 \N > 

= -i- — ° W gAG^eNr^xlsN . (7) 
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at tree- level in the Standard Model (Figure 2a). Here, J\ (./as) an d ^a (^as) denote the 
vector (axial vector) weak neutral currents of the quarks and electron, respectively 0. 
The anapole moment interaction of Eq. @) generates additional contributions to M^ V A 
when a photon is exchanged between the nucleon and the electron (Figure 2b). The 
corresponding amplitude is 



x 

Note that unlike iMy^, iM^ contains no (1 — 4sin 2 6 W ) suppression. Consequently, the 
relative importance of the anapole interaction is enhanced by 1/(1 — 4sin 2 @w*) r *~' 

10. This 

enhancement may be seen explicitly by converting Eqs. @ and into i?^ =0,1 : 

T=0 8V2na 1 a s 

lanapole" 1 - Asm 2 6 W g A { ) 



„t=1i 8V2na 1 a v 

lanapole q A 2 \- Asm 2 6 w g A K } 



The constants a SjV contain contributions from loops generated by the Lagrangians 
given in Section 3 and from counterterms in the tree- level effective Lagrangian of Eq. ([3j): 

a s ,v = a L s , v + a% ■ (H) 

In HBChPT, only the parts of the loop amplitudes non-analytic in quark masses can be 
unambigously indentified with a^ v . The remaining analytic terms are included in a° T v . In 
what follows, we compute explicityly the various loop contributions up through C(l/A 2 ), 
while in principle, a% should be determined from experiment. In Section 5, however, we 
discuss a model estimate for a° T v . 

Before proceeding with details of the calculation, it is useful to take note of the scales 
present in Eqs. @. The constants generally proportional to a product of strong 

and weak meson-baryon couplings. The former are generally of order unity, while the size 
of weak, PV couplings can be expressed in terms of g n = 3.8 x 1CT 8 , the scale of charged 



current contributions [IB|. One then expects the AM contributions to the axial radiative 



corrections to be of order 

Rr^- 'f:" A l T- a ^*-0.01 • (12) 

In some cases, the PV hadronic couplings may be an order of magnitude larger than g n . 
Alternatively, chiral singularities arising from loops may also enhance the AM effects over 
the scale in Eq. (|i~2"D. Thus, as we show below, the net effect of the AM is anticipated to 
be a 10-20 % contribution to R T r°'\ 
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III. NOTATIONS AND CONVENTIONS 



Since much of the formalism for HBChPT is standard, we relegate a detailed summary 
of our conventions to Appendix A. However, some discussion of the effective Lagrangians 
used in computing chiral loop contributions to a SjV is necessary here. Specifically, we 
require the parity-conserving (PC) and parity-violating (PV) Lagrangians involving pseu- 
doscalar meson, spin-1/2 and spin-3/2 baryon, and photon fields. For the moment, we 
restrict ourselves to SU(2) flavor symmetry and generalize to SU(3) later. The relativistic 
PC Lagrangian for tt, N, A, and 7 interactions needed here is 

C po = S-TrD'ED^ + N(iV^ - m N )N + g A NA^ l5 N 
+^N(c s + c v r 3 )a^F+N 

+y^LV75 + f (l»AV + A^ Ju ) l5 + | T ^7 a 757,]^ 
-ie^TfF+T» + [^-f£(d s + d v r 3 )rj 5 F+N + h.c] (13) 

X X 

where X> M is a chiral and electromagnetic (EM) covariant derivative, E = exp(zr • tt/F^) 
is the conventional non-linear representation of the pseudoscalar field, A is a nucleon 
isodoublet field, is the A field in the isospurion formalism, F^ v is the photon field 
strength tensor, and A^ is the axial field involving the pseudoscalars 

A, = + 0{^) (14) 

with being the EM covariant derivative. Explicit expressions for the fields and the 
transformation properties can be found in Appendix A. The constants c s , c v determined 
in terms of the nucleon isoscalar and isovector magnetic moments, ca is the A magnetic 
moment, d s , d v are the nucleon and delta transition magnetic moments, and Zq is the 
off-shell parameter which is not relevant in the present work Our convention for 75 
is that of Bjorken and Drell Ij| . 



In order to obtain proper chiral counting for the nucleon, we employ the conventional 
heavy baryon expansion of C PC , and in order to cosistently include the A we follow the 
small scale expansion proposed in [|17| . In this approach energy-momenta and the delta 
and nucleon mass difference S are both treated as 0(e) in chiral power counting. The 
leading order vertices in this framework can be obtained via P+FP+ where T is the original 
vertex in the relativistic Lagrangian and 
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are projection operators for the large, small components of the Dirac wavefunction respec- 
tively. Likewise, the 0(1/ m^) corrections are generally propotional to P + TP_/m^. In 
previous work the parity conserving irNA'y interaction Lagrangians have been obtained 
to 0(l/m 2 N ) |T7|]. We collect some of the relevant terms below: 



C p v c = N[iv ■ D + 2g A S ■ A]N - iTftiv ■ D 13 - 5 tJ 8 + gi 3 ■ A l3 ]T£ 



A 4 

1 ( ,-.. . 2 ,,2 



+- N{ (v ■ D) 2 — D + [S,, S„] [£>", D v \ 

2rriN 1 



-Wa 



(S ■ Dv A + v AS ■ D))N + ■■■ (16) 



where is the Pauli-Lubanski spin operator and 5 = m& — m^. 

The PV analog of Eq. (]T3|) can be constructed using the chiral fields X£ R defined 
in Appendix A and the spacetime transformation properties of the various fields in Eq. 
(|13|). We find it convenient to follow the convention in Ref. 0] and separate the PV 
Lagrangian into its various isospin components. The hadronic weak interaction has the 
form 

H w = ^J=JxJ X ] + h-c. , (17) 

where J\ denotes either a charged or neutral weak current built out of quarks. In the 
Standard Model, the strangeness conserving charged currents are pure isovector, whereas 
the neutral currents contain both isovector and isoscalar components. Consequently, TC W 
contains AT = 0, 1, 2 pieces and these channels must all be accounted for in any realistic 
hadronic effective theory. 

Again for simplicity, we restrict our attention first to the light quark SU(2) sector. (A 
general SU(3) PV meson-baryon Lagrangian is given in the Appendix and is considerably 
more complex.) We quote the relativistic Lagrangians, but employ the heavy baryon 
projections, as described above, in computing loops. It is straightforward to obtain the 
corresponding heavy baryon Lagrangians from those listed below, so we do not list the 
PV heavy baryon terms below. For the irN sector we have 

£at=o = h° v NA^N (18) 
£at =1 = ^-NYNT^A.Xl) - ^-NYl 5 NTr(A,X 3 _) (19) 



-.F^NX 6 N 



2^2 n 

£at= 2 = h 2 v l ab N[X a R A,X b R + X a L A,X b L ]YN (20) 
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The above Lagrangian was first given by Kaplan and Savage (KS) [iTj. However, the 
coefficients used in our work are slightly different from those of Ref. []TTJ] since our definition 
of An differs by an overall phase (see Appendix A). Moreover, the coefficient of the second 
term in the original PV AT = 2 NNtttt Lagrangian in Eq. (2.18) was misprinted in the 
work of KS, and should be 2h\ in their notation instead of hy as given in Eq. (2.18) of 

ED- 

The term proportional to h n contains no derivatives and, at leading-order in 1/F n , 
yields the PV NNtt Yukawa coupling traditionally used in meson-exchange models for the 
P V NN interaction [fT6|JT9[1 . The P V 7-decay of 18 F can be used to constrain the value of h n 
in a nuclear model-independent way as discussed in Ref. fll9l , resulting in h n = (0.7±2.2)g 7r 



15|j . Future PV experiments are planned using light nuclei to confirm the 18 F result. 



The coupling h v has also received considerable theoretical attention |lqJ2^J2Q|21[| and is 
particularly interesting since it receives no charged current contributions at leading order. 

Unlike the PV Yukawa interaction, the vector and axial vector terms in Eqs. (p!8|-^0|) 
contain derivative interactions. The terms containing h\, h\ start off with NNtttt inter- 
actions, while all the other terms start off as NNtt. Such derivative interactions have not 
been included in conventional analyses of nuclear and hadronic PV experiments. Conse- 
quently, the experimental constraints on the low-energy constants h v , h l A are unknown. 
The authors of Ref. []TTJ used simple dimensional arguments and factorization limits to es- 
timate their values, and we present additional phenomenological considerations in Section 
6 below. We emphasize, however, that the present lack of knowledge of these couplings 
introduces additional uncertainties into R^~ 0,1 . 

In addition to purely hadronic PV interactions, one may also write down PV EM 
interactions involving baryons and mesons^. The anapole interaction of Eq. ([|) represents 
one such interaction, arising at (9(1/A^) and involving no 7r's. There also exist terms at 
(9(1/A X ) which include at least one tt: 

C ^ PV = £L Na ^[ F ^ x!] + iV + ^-Na^F^N + ^AV^, Xl\ + N . (21) 

X XX 

The corresponding PV Lagrangians involving a N — > A transition are somewhat more 
complicated. The analogues of Eqs. (|18| - ^0D are 



£%=o = fie abc m l5 [X a L A,X b L + X a R A„X b R ]T? 

+ gi N[A„ X a _] + T£ + g 2 N[A„ X a _} + T» + h.c. (22) 

Cti=i = f2e ab *Ni lb [A^Xl] + T£ + f 3 e ab3 Ni l5 [A^X a + ]„T£ 

+g 3 N[(X a L A IM Xl - XlA.Xl) - {X%A ll X% - X 3 R A,X^ 
+g 4 {N[3XlA»(XlT^ + X 2 L T 2 ) + 3(^X1^ + X 2 L A»X 3 L T 2 ) 



^Note that the hadronic derivative interactions of Eqs. (p^-pO|) also contain 7 fields as required 
by gauge-invariance 
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-2(XlA»Xl + X\A*X\ - 2XlA^Xl)T^} - (L «-> R)} + h.c. (23) 

Cti=2 = Ue abd I cd Ni l5 [XtA^X b L + X R A„X b R ]T? 

+f 5 e ab3 Ni l5 {X a L A^Xl + XlA,X a L + (L <- R)]T» 

+g 5 l ab N[A,,X a _} + T!; + g 6 l ab NiA f ,,X«] + T» + h.c. , (24) 

where the terms containing /, and gi start off with single and two pion vertices, respec- 
tively. 

Finally, we consider PV 7AA interactions: 



C£ N = ie^f^F+N + ie^f^[F+ , X 3 + ] + N (25) 



A X A X 



A x A x 



+ze^f 3 V[F M - , Xi] + A + ie^f 3 ^[F- u , X!]_A + h.c. (26) 
A x A x 

The PV 7AA vertices di_ 3 , c? 4 _ 6 and dj^ 8 are associated at leading order in l/F n with 
zero, one and two pion vertices, respectively. All the vertices in (|l~8D -(p5|) are 0{p) or 
C(l/A x ) except h V) which is Yukawa interaction and of 0(p°). As we discuss in Appendix 
C, we do not require PV interactions involving two A fields. 



IV. CHIRAL LOOPS 

The contributions to a S)U arising from the Lagrangians of Eqs . (|I~8| - pO| ) are shown in 
Figure 3. We regulate the associated integrals using dimensional regularization (DR) and 
absorb the divergent — l/(d — 4) — terms into the counterterms, a° T v . The leading contribu- 
tions arise from the PV Yukawa coupling h n contained in the loops of 3a-f. To 0(1/A X ), 
the diagrams 3e,f containing a photon insertion on a nucleon line do not contribute. The 
reason is readily apparent from examination of the integral associated with the amplitude 
of Figure 3e: 

... . , V2g A f d D k i(S-k) 
%M 3e = xe^h^v ■ e — — 



(27r) D v ■ k v ■ (q + k) k 2 — m\ + ie 



< 2 ^9a c r° , f 1 , f d D k k 







where is the photon momentum, e is the photon polarization vector, s has the dimen- 
sions of mass, and we have Wick rotated to Euclidean momenta in the second line. From 
this form it is clear that iM 3e oc S ■ v = 0. The sum of the non- vanishing diagrams Figure 
3a-d yields a gauge invariant leading order result, which is purely isoscalar: 
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a L a (Yl) = • (28) 

24 m n 

As the PV Yukawa interaction is of order 0(p°), we need to consider higher order 
corrections involving this interaction, which arise from the 1/mjv expansion of the nucleon 
propagator and various vertices. Since P + ■ 1 • P_ =0, there is no 1/mjy correction to the 
PV Yukawa vertex. From the 1/m.jv iVA^ terms in Eq. JIB]) we have 



<(Y2) = ^e^Aln(^) 2 , (29) 

where [i is the subtraction scale introduced by DR. Finally, the 1/rriN correction to the 
strong ttNN vertex, contained in the term oc g& in Eq. flT3|), yields 



<(F3) = -^.egA^H—) 2 ■ (30) 
487T m N m n 



These terms are also isoscalar, and the results in Eqs. ( p8| - |30| ) are fully contained in the 
previous analyses of Refs. |5|,[7|,P1| . 



For the interactions in Eqs. (IIBH20) containing hy, the eight diagrams Figure 3a-h 



must be considered. Their contribution is purely isovector — 

<{V) = \eg A {h v + \h 2 v ) ln(^) 2 . (31) 

O 777/ 

— and was not included in previous analyses. 

The contribution generated from the two-pion PV axial vertices in Eqs. (|19l - p0|) comes 
only from the loop Figure 3i and contains both isovector and isoscalar components: 

<{A) + a L v (A)r 3 = -\e{h\ + h\r 3 ) ln(^) 2 . (32) 

O 777 Tj- 



a result first computed in Ref. |l I 



In principle, a variety of additional contributions will arise at 0(1/K^). For exam- 
ple, insertion of the nucleon magnetic moments {i.e. the terms in Eq. (^) containing 
c s>v ) into the loops Figure 3e,f — resulting in the loops of Figure 5a,b — would in principle 
generate terms of (^(l/A 2 ,) when the PV Yukawa interaction is considered. As shown 
in Appendix C, however, such contributions vanish at this order. Similarly, the entire 
set of A intermediate state contributions shown in Figure 4, as well as those generated 
by C],y and £],y N in Figure 6, vanish up to 0(1/ A 2 ). The reasons for the vanishing of 
these various possible contributions is discussed in Appendix C. Thus, the complete set 



of SU(2) loop contributions up to 0(1/A^) are given in Eqs. ([28|]3 



Because m c — m s » m s — m u ^ and A x >> m s , it may be appropriate to treat the 
lightest strange and non-strange hadrons on a similar footing and extend the foregoing 
discussion to SU(3) chiral symmetry. A similar philosophy has been adopted by several 
authors in studying the axial charges and magnetic moments of the lightest baryons 
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I4J22H25] . In what follows, we consider the possibility that kaon loop contributions, 
introduced by the consideration of SU(3) symmetry, may further enhance the anapole 
contribution to Ra- 

Before proceeding along these lines, however, one must raise an important caveat. 
When kaon loop corrections are included in a HBChPT analysis, higher order chiral 
corrections may go as powers of m K /A x ~ 0.5. Consequently, the convergence of the 
SU(3) chiral expansion remains a subject of debate [p6fl . Fortunately, no such factors 



appear in the present analysis through 0(1/A X ) so that at this order, we find that kaon 
loop effects in Ra are generally tiny compared to those involving pion loops. Whether or 
not higher-order terms (e.g., those of 0(1/ A x x m K /A x ) contribute as strongly as those 
considered here remains a separate, open question. 

To set our notation, we give the leading strong-interaction SU(3) Lagrangian. Since 
the K° and 77 are neutral, loops containing these mesons do not contribute to the AM 
through 0(1/ A x ) and we do not include their strong couplings below. For the proton 
the possible intermediate states are Y,°K + , AK + while for the neutron only T,~K + can 
appear. The necessary vertices derive from 

C = 2g A NS ■ AN + 2g NAK [(NS ■ K)A + A(S ■ K ] N)] (33) 
+2g NSK [S ■ K^N + NZS ■ K] , 

where g NAK = -[(1 + 2a)/V&\g A , 9nzk = (1 - 2a)g A with g A = D + F, a = F/(D + F) 
and D, F are the usual SU(3) symmetric and antisymmetric coupling constants. 

The general pesudoscalar octet and baryon octet PV Lagrangians are given in the 
Appendix B. They contain four independent PV Yukawa couplings, 20 axial vector cou- 
plings (/i^-type) , and 22 vector couplings (/iy-type). For simplicity, we combine the SU(3) 
couplings into combinations specific to various hadrons — e.g. the leading PV Yukawa in- 
teractions are 

^Yukawa = ~ ih ^(P n ^ + ~ npir~) - ih pTi o K (pY? K + - t°pK~) 
-ih nE - K (nJ:-K + - YrnK") - ih pAK (pAK + - ApK~) + ■■■ . (34) 

In terms of the SU(3) couplings listed in Appendix B, the Hbbm have the form 

K = -2V2(hx + h 2 ) 
h p zo K = -[/ii - h 2 + ^(^3 - h 4 )] 
h n T.-K = ^2h V Y$ K 

' • (35) 



h 



pAK 



+ V3h 2 + h 3 + 3h 4 

v3 



Similarly, we write for the vector PV interaction 

h?F- K+ h pAK+ 



^—n^H-D^K + - -t—p^AD^K + + h.c. + • • • , (36) 



10 



and for the axial PV two pion and kaon interactions 

, pir ,pK 

7T 7T 
L. T17T LTlK 

+i-^nYl 5 n^ + D^- - tt-D^) + i-^n^ lb n{K + D ^ - K~D^K+) 

"it "it 

+ •••• (37) 

Expressions for these PV vector and axial coupling constants in terms of SU(3) constants 
appear in Appendix B. For illustrative purposes, it is useful to express the nucleon-pion 



couplings in terms of the h l VA of Eqs. (p^-pO|) for the SU(2) sector 



h p r + = k + \h\ 
h^ = h\+h\ 

h7 = h\-h\ . (38) 

The leading order contributions to a SjV arise only from the loops in Figure 3 where 
a photon couples to a charged meson. The charged kaon loop contributions to the a SjV 
can be obtained from the corresponding formulae for the 7r-loop terms by making simple 
replacements of couplings and masses. For example, for the PV Yukawa interactions, 
these replacements are: (a) for the proton case, m n — > m K , h n — > h p ^o K , g A 9nt,°k+ — 
9n?:k/V2 for TPK + intermediate states and h n — > h p \K, 9a — >■ 9nkk for AK + ; (b) for the 
neutron case, h n — > h n %-K, 9a —* 9ny.-k+ = 9ny,k for T,~K + intermediate state for the 
neutron case. Similar replacements hold for the vector PV coupling contributions. For the 
axial PV two-pion contribution we need only make the replacement h A n — > h A , — > itlk- 

Upon making these substitutions, we obtain the complete heavy baryon loop contri- 
bution to <D(l/A 2 x ) in SU(3): 

a L s = §9AhA-^ + -^H^) 2 ] 
24 m n 7r rrijsi m v 

v3 ,„ t n N , r A x 3 A x 1 / [i 



1 + 2a)g A h P AK[ + ln(- 



144 ttik is tuk 

+^§(l-2a)g A h nS - K [-^ + -^H >' 



32 rrix tcttin m K 

~h h A + h?) H~) 2 - Uhf + Nf) ln(^) 2 
o m n o rrix 



~(1 + 2a)g A h^ \n(-^r (39) 
72 rrix 



< = -7^r( 1 - 2 <*)9aKe-k[ — - + — - H 



96 m K Ttm N ttik' 
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73 
144 



(1 + 2a)g A h 



A v 3 A, 



pAK[ 



+ 



ln( 



rriK Trmjs! tjik 
~\{h P 2 - h n /) ln(^) 2 - Uhf - h A K ) ln(^) 2 

IB?,- O 777^ 

1 

6 i 



-lg A hr + H — ) 2 



-i t,P T -° K+ 

-^(l + 2a)^r + ln(^) 2 . 
72 rriK 



ln( 



A* x2 



(40) 



V. LOW-ENERGY CONSTANTS AND VECTOR MESONS 

A pure ChPT treatment of the anapole contributions to R A would use a measurment 
of the axial term in A LR (ep) and A LR (QE), together with the non-analytic, long-distance 
loop contributions, a^ v , to determine the low-energy constants, a°*. In the present case, 
however, we wish to determine whether there exist reasonable hadronic mechanisms which 
can enhance the low-energy constants to the level suggested by the SAMPLE results. 
Thus, we attempt to estimate a° T v theoretically 

Because they are governed in part by the short- distance (r > 1/A X ) strong interaction, 
are difficult to compute from first principles in QCD. Nevertheless, experience with 
ChPT in the pseudscalar meson sector and with the phenomenology of nucleon EM form 
factors suggests a reasonable model approach. It is well known, for example, that in the 
0(p A ) chiral Lagrangian describing pseudoscalar interactions, the low-energy constants 
are well-described by the exchange of heavy mesons p7| |. In particular, the charge radius 



of the pion receives roughly a 7% long-distance loop contribution, while the remaining 
93% is saturated by t-channel exchange of the p°. Similarly, in the baryon sector, disper- 
sion relation analyses of the isovector and isoscalar nucleon electromagnetic form factors 
indicate important contributions from the lightest vector mesons ||28|| . Thus, it seems 
reasonable to assume that t-channel exchange of vector mesons also plays an important 
role in the short- distance physics associated with the anapole moment. 

With these observations in mind, we estimate the coefficients af T v in the approximation 
that they are saturated by t-channel exchange of the lightest vector mesons, as shown 
in Figure 7. Here parity-violation enters through the vector meson-nucleon interaction 
vertices. We also use a similar picture for the electromagnetic nucleon form factors to 
determine the overall phase of a c s T v in the vector meson dominance approximation. To 
that end we require the PC and PV vector meson-nucleon Lagrangians |L6| : 



C% N = 9 p nnN[i, + «p^-]t " P"N (41) 
C P U C NN = gu, NN N[% + kJ^-]u^N (42) 
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'PC _ at r | i&fii/Q 



^nn = + k^-^-^N (43) 

and 

= N^PIK + (h° p + i)r 3 ]iV (44) 

CZnn = N-fwM + hlr 3 ]N (45) 
£+nn = NYlsMK + fy*)N . (46) 

(Note that we have adopted a different convention for 75 than used in Ref. ||16|| .) The 
coupling constants h p ^^ were estimated in Refs. [[T^,[^J and have also been constrained 
by a variety of hadronic and nuclear parity- violating experiments (for a review, see Ref. 

0)- 

For the V — 7 transition amplitude, we use 

= ^-F^ , (47) 

*JV 

where e is the charge unit, fy is the j-V conversion constant (V = p°, u>, 0), and V pv is the 
corresponding vector meson field tensor. (This gauge- invariant Lagrangian ensures that 
the diagrams of Figure 7 do not contribute to the charge of the nucleon.) The amplitude 
of Figure 7 then becomes 

h 1 A h° A h° A 

a c s T (VMD) = -/(^) 2 + ^(^) 2 + -^) 2 , (48) 
Jp rn p f u m u m 



«TW = ; '° + ? /V ~\ ±f + + . (49) 

Jp rn p f u m u 

The parity violating rho-pole contribution was first derived in [|5],[7| . However, the relative 
sign between h pN and f p is undetermined from the diagram of Figure 7 alone. Never- 
theless, we can fix the overall phase using two phenomenological inputs. Parity violating 
experiments in the p-p system constrain the sign of the combination g p is[h pN [j3l]|, |3"l"| , |H| | 



In particular, the scale of the longitudinal analyzing power, Al, is set by the combination 
of constants 

A L oc g pNN {2 + k v ) [h° p + h\ + h 2 p /V6] + g„ NN (2 + k s ) [h° u + h x u \ , (50) 

where the constant of proportionality is positive, Ky = 3.7 and Ks — —0.12. Using 
the standard values for the strong VNN couplings, one finds that Al has roughly the 
same sensitivity to each of the h l v appearing in Eq. ([50]) (modulo the 1 / y/6 coefficient of 
h 2 p ). From the 45 MeV experiment performed at SIN ||32|| , for example, one obtains the 



approximate constraint 15 



h° p + h l p + h 2 p /VQ + tii + hl~-2%±A , (51) 
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where the hy have are expressed in units of g n and where a positive sign has been assumed 
for Qvnn- Given this constraint, it is very unlikely that the product (h° p + h 2 /y/6)g P NN > 
unless the corresponding products involving h p and hPj 1 in Eq. (|50"D obtain anomalously 
large, negative values. In fact, a fit to hadronic and nuclear PV observables in Ref. [[0"J 
strongly favors a phase difference between the strong and weak VNN couplings. 
Experimentally, one also knows the isovector nucleon charge radius 



r 2 )^ = 6^|^|^o>0, (52) 



where 



<M = \0)\p>=eu(p')[F 1 (q 2 ) + l -^F 2 ( q 2 )}u(p) . (53) 

One may reasonably approximate the p° contribution to (r 2 ) T=1 using VMD ||28|| . 
The calculation is the same as above but with the weak hadronic coupling replaced by 
the strong coupling. The result is 

F f { f ) = il^-f— . (54) 
f p q 2 -m 2 p 

Then we have 

dF ™°tf) \ 9 P NN _ (55) 

dq 2 q ° f p m 2 p ' 

Comparing Eqs. ( [5*2] ) and ([55]), and noting that the p° generates a positive contribution 
to (r 2 ) T=1 [f28f| , we arrive at g P NN/ fp < 0. Combining this result with g P NNh p < as 
favored by the pp experiments |3T||3T] , |i"9fl we obtain the relative sign between h l p and f p : 
h„/ f p > 0. Accordingly we determine the relative signs for PV u, 0-nucleon coupling 
constants. 



VI. THE SCALE OF R A 



Expressions for the anapole contributions to and 1 in terms of the a StV 
5[). We may now use these expressions, along with the results in Eqs. 



=0,1 1 



appear in Eq. (["]). 

( |39|]4"0| ) and ( "4i^ - |49|) , to obtain a numerical estimate for the R^ 
use the global fit value for the weak mixing angle in the on-shell scheme, sin 2 6 



To do so, we 



g A = 1.267 ±0.004 
= A, 



5.26 @, U = 17, 



13 34 



a 



x . We express all the PV coupling constants in units of g n 
traditionally done [E~|,|n"j]. We obtain 



anapole - 

w = 0.2230 
F/(D + F) = 0.36, 
= 3.8 x 



10~ 8 as is 
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^ = °lanapole = ^~ 2 {0-^K + h\- 0.0036/i nS - K 
-0.033(/# 



+ -^=H + 0.2(h p A K + Kf ) - 0.006/i P AAr 



V2 

+0.088h p v AK+ - 0.26\h\\ - 0.08|/i°| - 0.05|/i°|} 



(56) 



4. 



^lanapole = ~ 0.6(4 + ^K) ~ 0.0012/i nS - K - 0.033(-/i 



V 



v/2 



+ 0.2(/i^f - - 0.006/w^ + 0.088^ AK+ 



-0.26(|^| + 



\h 2 \ 



0.087|/£| + 0.05|M|} 



(57) 



where we have set the phase of the vector meson contributions as discussed above, and 
used the relations in Eq. (|38|). 
The expressions in Eqs. 



57 ) illustrate the sensitivity of the radiative corrections 



to the various PV hadronic couplings. As expected on general grounds, the overall scale 
of R^ =0,1 is at about the one percent level [see Eq. (0)]. In terms of the conventional PV 
couplings, R^ =0 is most sensitive to h n and h 1 , while R^ =1 is most strongly influenced by 
h° p + h 2 /y/6. The corrections also display strong dependences on the couplings h l VA not 
included in the standard analysis of nuclear and hadronic PV. In particular, the couplings 
h\ and h° v + Ah 2 v /3 are weighted heavily in 1 . In general, the sensitivity to the 
PV NYK couplings is considerably weaker than the sensitivity to the NNtt and NNp 
couplings. 

In order to make an estimate of R^ =0,1 , we require inputs for the PV couplings. To 
that end, we use the "best values" for h n , h p , and h l u given in Ref. p9| . These values 
are consistent with the fit of Ref. [||J. For the hi 
The analyses given in Refs. 



16,19,29 



we use the "best values" of Ref. fTB]] . 
together with experimental input, also allow for 
the standard couplings to take on a range of values. For example, the ranges for the 
given in Refs. |l6lp9j correspond to 



- 33 < h°, + h}„ < 13 



(58) 



In order to maintain consistency with the experimental constraint of Eq. (^T|), one then 
requires 



< h° + h\ + h 2 p /V6 < -45 



(59) 



29 . Indeed, 



We adopt this range even though it is smaller than the range given in Ref. 
allowing the h p to assume the full ranges given in Ref. [2~9~| would require the h 1 ^ to vary 
outside their corresponding theoretical "reasonable ranges" if the constraint of Eq. (|5T|) 
is to be satisfied. Since one expects << \h° p ,2 \ [fT6| , p9|1 , we have a reasonable range 
of values for the important isoscalar p contribution in Eq. ([57]), and the rather broad 
range of values allowed for the h % contributes significantly to our estimated uncertainty 
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in -R^ -1 . For h^ ^, we use the ranges of Refs. |1|^9]]0. In contrast to the situation with 



the h? contribution, however, the variation in the ^ over their "reasonable ranges" has 
negligible impact on our estimated theoretical uncertainty. 

Estimating values for the Yukawa couplings h^YK and for the hy t A is more 



problematic — to date, no calculation on the level of Ref. |16| has been performed for 



such couplings. Estimates for hy based on dimensional and factorization arguments, 
were given in Ref. [|ll[] and generally yielded values for hy,A in the non-strange sector on 
the order of g n . For our central values, then, we take hy A = g n , resulting in roughly 
1% contributions from the PV vector and axial vector interactions. Without performing 
a detailed calculation as in Ref. [16], one might also attempt to determine reasonable 
ranges for these parameters by looking to phenomenology. To that end, the authors of 
Ref. fll| considered analogies between the axial vector PV operators of Eq. (|l9|-p0|) and 



contact operators needed to explain the size of AI = 1/2 hyperon P-wave decay ampli- 
tudes. From this analogy, these authors conclude that \h A \ ~ lOg^ may be reasonable. 
However, whether such large ranges are consistent with nuclear PV data remains to be de- 
termined. In the absence of such an analysis, which goes beyond the scope of the present 
work, we adopt the range —I0g n < h\ < I0g n suggested in Ref. (JTI]]. The corresponding 
uncertainties in the -R^ -0 ' 1 are roughly ±10%. 

The implications of phenomenology for the hy are even less clear than for the h\. 
However, we note that large values hy ~ ilOgv do not appear to be ruled out by hadronic 
and nuclear PV data. At tree-level, for example, the vector terms in £at=o i 2 do not 
contribute to the PV NN interaction through the one 7r-exchange amplitudes of Figure 
8a. It is straightforward to show that the corresponding amplitude vanishes for on-shell 
nucleons |. Thus, at this level, purely hadronic PV processes are insensitive to the h v 
and provide no constraints on these couplings. In PV electromagnetic processes, however, 
the h v do contribute through PV two-body currents, such as those shown in Figure 8b. 
Nevertheless, one expects the impact of PV two-body currents to be considerably weaker 
than that of the PV NN potential. The PV 7-decay of 18 F, for example, is dominated by 
the mixing of a nearly-degenerate pair of (J 7r , T) = (0~,0) and (0 + , 1) states. The small 
energy denominator associated with this parity-mixing enhances the relative importance 
of the PV NN potential by roughly two orders of magnitude over the generic situation with 
typical nuclear level spacings. By contrast, the PV two-body currents do not participate 
in parity-mixing and receive no such enhancements. A similar situation holds for PV 
electromagnetic processes in other nuclei of interest. Hence, we expect the PV 7-decays 
of light nuclei to be relatively insensitive to the h l v , even if the latter are on the order of 
10g n . Consequently, we rather generously take —lOg-^ < h v + 4/iy/3 < 10g n , yielding a 



* Allowing the h l u to assume positive values would require a sign change on the correspnding 
terms in Eqs. ( |56|J)7| ). 

§The on-shell approximation is generally used in deriving the PV NN potential from Feynman 
diagrams. 
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±7% contribution to the uncertainty in R^ =1 - Allowing similarly large ranges for the PV 
NYK couplings has a negligible impact on the uncertainty in the R^~ 0,1 . 

With these input values for the PV couplings, we arrive at the anapole contributions 
to R^ =0,1 shown in Table I. The latter must be added to the one-quark Standard Model 
contributions, also shown in Table I. We compute the one-quark corrections using the 
on-shell parameters given in Refs. P,|5B|. We emphasize that the quoted values for the 
Ra~ 0,1 are renormalization scheme-dependent. The relative size of the isovector one-quark 
corrections are smaller, for example, in the MS scheme, where one has i?J =1 (SM) = 
—0.18 and i?^ =0 (SM) = 0.07. The corresponding tree- level amplitude, however, is also 
smaller by a factor of ~ 1.44 than the on-shell tree- level amplitude. A reader working 
in the MS scheme should, therefore, take care to adjust the tree-level amplitude and 
SM radiative corrections appropriately from the on-shell values used here. Moreover, the 

(T) 

anapole contributions to the R A will be a factor of 1.44 larger in the MS scheme since 
the tree-level amplitude is correspondingly smaller]**]. 

Adding the one-quark and anapole contributions yields a large, negative value for 
i?J =1 . This result contains considerable theoretical uncertainty, mostly due to our liberal 
assignment of reasonable ranges to the hy,A- Even with this generous theoretical uncer- 
tainty, however, R^ =1 is still roughly a factor of two away from the apparent SAMPLE 
result. Compared with the one-quark SM contribution, the many-quark anapole contribu- 
tion is relatively small - though it does push the total in the right direction. The isoscalar 
correction, -Rj =0 , is considerably smaller in magnitude than R^ =1 yet retains a sizeable 
theoretical uncertainty. 



**Note that gives the ratio of the isoscalar, axial vector amplitude to the tree-level 

isovector, axial vector amplitude. The sign of R^ =0 as defined here is opposite that of Ref. ||]. 
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TABLES 



Source 


R T A =1 


R T=0 


One-quark (SM) 


-0.35 


0.05 


Anapole 


-0.06 ±0.24 


0.01 ±0.14 


Total 


-0.41 ±0.24 


0.06 ± 0.14 



TABLE I. One-quark Standard Model (SM) and many-quark anapole contributions to 
V{A) x A(N) radiative corrections. Values are computed in the on-shell scheme using 
sin 2 9 W = 0.2230 . 
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VII. CONCLUSIONS 



In view of the preliminary SAMPLE result for PV quasielastic electron scattering from 
2 H, we have up-dated our previous calculation of the axial vector radiative corrections 
i?^ -0,1 . Using the framework of HBChPT, we have computed all many-quark anapole 
contributions through 0(1/A 2 ). We include new one- loop contributions involving the PV 
vector couplings, h % v and estimate the scale of the analytic, low-energy constants using 
resonance saturation. We fix the sign of the latter using the phenomenology of PV pp 
scattering and of nucleon EM form factors. We also show that large classes of loops 
involving decuplet intermediate states, magnetic insertions, and PV EM insertions vanish 
through 0(1/A 2 ). Finally, we extend the previous analyses to include SU(3) symmetry, 
and determine that the impact of kaon loops is generally negligible. In the end, we 
find that R T A =1 — though large and negative — is still a factor of two or so away from the 
suggestion that R^ =1 ~ — 1 from the SAMPLE experiment. Even allowing for considerable 
theoretical uncertainty - dominated by the PV couplings h l VA - there remains a sizable 
gap between our result and the preliminary experimental value. 

There exist a number of possible additional contributions to i?^ =0,1 not considered 
here which may ultimately account for the apparent experimental result. The most ob- 
vious include higher-order chiral corrections. This appears, however, to be an unlikely 
source of large contributions. On general grounds, we expect the size of the 0(1/A^) con- 
tributions to be suppressed by m/A x relative to those considered here, where m denotes 
a pseudoscalar mass. For kaon loops, this suppression factor is only ~ 1/2; however, at 
(9(1/A 2 ) kaon loops generate at most a few percent contribution to i?^ =0,1 . The suppres- 
sion factor for the next order pionic contributions is closer to 0.1. Hence, it would be 
surprising if the next order in the chiral expansion could close the factor of two gap with 
experiment. 

More promising sources of sizeable contributions include Z — 7 box graph contribu- 
tions, where the full tower of hadronic intermediate states is included, as well as parity- 
mixing in the deuteron wavefunction. At a more speculative level, one might also consider 
contributions from physics beyond the Standard Model. For example, the presence of an 
additional, relatively light neutral gauge boson might modify the SM V(e) x A(q) ampli- 
tudes and contribute to R T A =1 . A popular class of Z' models are generated by E 6 symmetry 



3q| . The contribution of an extra, neutral weak ~Eq gauge boson Z' is given by 



^I =1 (new) = | 1 sin 2 0^ , (60) 
5 1 — 4sm 6 W G M 

where is a mixing angle which governs the structure of an additional U(l) group in Eq 
theories |I6| and G'^ is the Fermi constant associated with the new U(l) group |37|. Note 
that this contribution has the wrong sign to account for the large negative value of R^ =l - 
Alternatively, one might consider new tree-level interactions generated by supersym- 
metric extensions of the SM. Such interactions arise when R-parity, or equivalently, B — L, 
is not conserved (B and L denote baryon and lepton number, respectively). The contri- 
bution from R-parity violating SUSY interactions is given by |38| , |37| , |39 
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.R^new) 



where x = sin 2 9 W , 



' 1 A' llfc (4)-A' ljl (gl)-A m (el)(l-4x + 4A x )l , (61) 



I- Ax 



x(l — x) ( 1 \ , . 

A* = — r- ~0.3 , 62 

Ar is a radiative correction, and where 

^ = 4755$? ' (63) 

with / denoting the superpartner of fermion / and i, j, A; labeling fermion generations. 
The terms having a prime are semileptonic whereas the un-primed terms are purely lep- 
tonic. In principle, the correction in Eq. ( |5"T| ) could generate a negative contribution to 
R^ =1 - However, the various other electroweak data constrain the terms appearing in this 
expression. For example, relations between and other SM parameters require 

- 0.0023 < A m (4) < 0.0028 , (64) 

at 90 % C.L., so that the first term in Eq. (|6l| ) cannot provide the large negative contri- 
bution needed to explain the SAMPLE result. Similarly, assuming only the semileptonic 
R-parity violating interactions modify the weak charge of nuclei, the recent determination 
of the cesium weak charge by the Boulder group PD| , |TU| implies that 



0.0026 < 2.6A' m (4) - 2.9A' lj . 1 (g{) < 0.015 , (65) 

at 95 % CL. (for m# = 300 GeV). Thus, it appears unlikely that the second term in Eq. 
( pip could enhance R^ =1 by a factor of two. 

In short, two of the most popular new physics scenarios having implications for low- 
energy phenomenology appear unlikely to enhance R^ =1 significantly. Thus, if more 
conventional hadronic and nuclear processes cannot account for the SAMPLE result, one 
may be forced to consider more exotic alternatives. 
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APPENDIX A: FORMALISM 



In this section we first review the general parity and CP conserving Lagrangians 
including N, 7T, A, 7 in the relativistic form. We follow standard conventions and introduce 

£ = e 2 , Z = , vr = \ a r a (Al) 

with F,,- = 92.4 MeV being the pion decay constant. The chiral vector and axial vector 
currents are given by 

a, = -k^e f - e f ^o = -^f + o(vr 3 ) 

Z ±< n 
V» = \(£D^ + ^D,0 . (A2) 

and we require also the gauge and chiral covariant derivativs 

Df/ir = d^TT - ieA^[Q, 7r] 
V li = D lt + V fl , (A3) 



with 



«=(!-j) < A4 > 



and ./L being the photon field. The chiral field strength tensors are 



with 



<? = (Jo) < A6 > 

acting in the space of baryon isodoublets. 

For the moment, we restrict our attention to SU(2) flavor space and consider just it, 
N, and A degrees of freedom. We represent the nucleon as a two component isodoublet 
field, while for the A, we use the isospurion formalism, treating the A field T^(x) as a 
vector spinor in both spin and isospin space flT7] with the constraint r*T^(a;) = 0. The 
components of this field are 



A+\ m+ /A++ \ /A 



^-fi[l-l •V'&wl = ' (A7) 

The field T l also satisfies the constraints for the ordinary Schwinger-Rarita spin-| field, 
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7 M T; = and p% = . (A8) 

We eventually convert to the heavy baryon expansion, in which case the latter constraint 
becomes v = with the heavy baryon velocity. 

It is useful to review the spacetime and chiral transformation properties of the various 
fields. Under a chiral transformation, 

X? M -> C/P^C/t, (A9) 

and 

N —> UN , — ► UTp , £ -f L£i? f etc. (A10) 

In the 577(2) sector parity violating effects are conveniently described by introducing 
the operators fill : 

X£ = £V£, AS = ^e t , A^ = A2±A^ . (All) 

which transform as 

xU^uxhtrt . (A12) 

with the index a rotating like a vector of SU(2)l and SU{2)r respectively. 

The P and CP transformation properties of these fields are shown in Table 2. 



Field 


p 


CP 


\ 
N 

xi 

F ± 


7o iV 


7oCA> T 

5{a)Xl a 
8{a)Xp 



TABLE II. Parity (P) and CP transformation properties of chiral fields. Here, T denotes 
the transpose, C is the charge conjugation matrix (C = ^7270 in the Dirac representation) and 
S(i) = l,i = 1,3 and 8(2) = -1. 



Finally, we note that in the Lagrangians of Section III, one has the following defini- 
tions: 

T% = $*Vp - 2ie ijk V k 
«J = Tr[r%] 

= ^3/2^M^3/2 ' (^13) 

where ^ 2 = — \^ijhT k is the isospin 3/2 projection operator. 
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APPENDIX B: THE SU(3) PARITY VIOLATING AND CP CONSERVING 

LAGRANGIAN 



In this Appendix we list the parity violating and CP conserving 577(3) Lagrangian for 
the pseudosclar meson octet and baryon octet. We are interested in the diagonal case of 
the parity violating electron nucleon scattering. Hence, we include only those interaction 
terms that ensure strangeness and charge conservation at each vertex. In the following we 



use f = e&, 7i = ±7r«A a , X£ = X% = X a ± = X a L ±X a R , [A, B} ± = AB ± BA. 



We classify the parity violating Lagrangian according to isospin violation AT = 0, 1, 2, 
which arises from the operators of X£, X^, their combinations and products. The AT = 2 
piece comes from the operators X ab {X10X b L ± (L <-> R)} with O — N, N, A^ and 



where a, b = 1,2,3. Several operators contribute to the AT = 1 part, like 
Xl,f ab {X a L OX b L ± (L <-> R)},d 3ab {X£OX b L ± (L <-> R)} where f ab ^d abc are the an- 
tisymmetric and symmetric structure constants of SU (3) algebra. With the requirement 
that the final Lagrangian be hermitian, parity-violating and CP-conserving, the operator 
with f 3ab vanishs. For the AT = part relevant operators are 1, X|, f 8ab {XlOX b L ±(L <-»■ 
R)},d 8ab {X£OX b L ± (L <-> R)} } 5 ab {X a L OX b L ± (L <-> i?)}. For the same reason the / 8ab 
structure does not contribute. Note the matrix identity \ a X b \ a = 4(C 2 (3) — |C2(8))A 6 , 
where 6*2(3), 6*2(8) are the Casimir invariants of the basic and adjoint representations 
of 577(3) group respectively. Hence, the operator containing 5 ab is identical to the unit 
operator. 

Based on these considerations, we obtain 




(Bl) 



£Xt=o = h 3 F w TrN[Xi, N} + + h 4 F n TrN[X 8 , 7V]_ + v^rN^A^ N\. 

+v 2 TrN^[A lx ,N}_ + — TrN^A^NXl + —TrN^X 8 + NA i 



+ j. T rNY[Xl, A,] + N + -fTrN^N[X 8 + , A,} + + a.TrN^A.NX 8 
+a 6 TrN 1 » l5 X 8 _NA fl + a 7 TrNYl 5 [X 8 _, AJ+N + a 8 TrN^ l5 N[X 8 , A^} + 



+V3v u d 8ab Tr{NYNX a L A^X b L + {L ^ R)} 



+V3v 12 d 8ab Tr{NYX a L A fl X b L N + (L <-> R)} 
+V3v l3 d 8ab Tr{NYX a L N[X b L , + (L <-> R)} 
+VZv 14 d 8ab Tr{NY[X a L , A^] + NX b L + (L <-> R)} 
+V3a 9 d 8ab Tr{NYl 5 NX a L A fl X b L - (L <-> R)} 
+V3a l0 d 8ab Tr{NYl 5 X a L A^X b L N - (L <- i?)} 
+v / 3a 11 rf 8a6 Tr{iV7^7 5 X2iV[X^,^] + - (L <- i?)} 
+v / 3ai 2 rf 8afc Tr{iV7^7 5 [X£,^] + iVXi - (L <- R)} , 



(B2) 
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PV 



/iii^TriVfXi, N] + + h 2 F n TrN[X 3 , N}_ + -fTrN^A^NX; 



AT=1 — 



+ 



+ -t T rNrX 3 + NA, + -jTrNY[X 3 + , AJ+N + -^TtN^N[X 3 + , A,] + 

+a l TrNYl 5 A fl NX 3 _ + a 2 TrN^ l5 X 3 _N A^ + a 3 TrNYl 5 [X 3 , A^+N 
+a 4 TrN^ 5 N[X 3 _, A^] + + v^d 3ab Tr{NYNX a L A^X b L + (L <-> R)} 

+v 1& d 3ab Tr{N^X a L A^X b L N + (L <-> R)} 
+v l7 d 3ab Tr{NYX a L N[X b L , A^ + (L <- R)} 
+v 18 d 3ab Tr{NY[X a L , A^] + NX b L + (L <- R)} 
+a 13 d 3ab Tr{NYl 5 NX a L A^X b L - (L R)} 
+a 14 d 3ab Tr{NYl 5 X a L A ll X b L N - (L <-> R)} 
+a 15 d 3ab Tr{Nri 5 X a L N[XlA II } + - (L <- R)} 
+a 16 d 3ab Tr{NYlrAX a L , A,] + NX b L - (L <- R)} , (B3) 

£l¥=2 = V -fT ab Tr{N^NX a L A^X b L + (L <-> R)} 
+ —X ab Tr{NYXlA^X b L N + (L <-> R)} 
+^T*Tr{N-fXlN[XlA ll ] + + (L - i?)} 

+ °fX^Tr{N^ lb XlN[XlA,\ + - (L <-> R)} 

+^ T *Tr{N'fT B [Xl, A,} + NX b L - (L ~ R)} . (B4) 



These Lagrangians contain 4 Yukawa couplings, 20 axial vector couplings and 22 
vector couplings, all of which should be fixed from the experimental data or from model 
calculations. In reality, however, we have only limited information which constrains a few 
of them. It is useful to expand the above Lagrangians to the order involving the minimum 
number of Goldstone bosons and to collect those vertices needed in the calculation of Ra- 



C-Yukawa = 2"\/2i(/ii + h 2 ){pmr + - fijm~) 
+t[h l -h 2 + V3{h 3 - h A )]{pY?K + - t° P K~) 
+V2i[h -h 2 + V3(h 3 - h A )}(nT 1 -K + - 



+ V3h 2 + h 3 + 3h 4 ](pAK + - KpK~) + ■■■ . 



(B5) 



fly 



F 



•+ 




F 



■p^°D^K + 
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v -—n 1 *YTD ll K + - -V— Pl »AD lx K + + h.c. + ■■■ , (B6) 

-f-rv t 'it 



where 



hT + = ^ + ^f(vi4 ~ V 12 ) + f(v 7 + V 9 ) + ^V 20 

hf° K+ = l( Vl -v 2 + v 4 + v 6 ) + v -^f + |(un - via - vis - «2i - \vn) + 2v l8 
hf-~ K+ = j-(vi -v 2 + v 6 - v 4 ) + ±(v s - Uio) + ^(Ul7 + v 2 i) + ^(uii - via - v 15 ) 

, P AK+ _ 1 / . 2,, 4„ , 16,, 2„ 2,, , 4„ 17,, 

fty — ^( 2" + 3 Vn ~ 3 Vl2 + Y Vl3 ~ 3 Vl4 ~ 3 Vl5 + 3^ 16 _ T^ 17 

+ iv 18 -lv 2 + f-v 5 + f) + ^^ + ^^ . (B7) 



/i p7r 

J 7T 

UPK 

J TV 

UnK 

-i^nYlMK+D^K- - K-D^K + ) + • • • , (B8) 

J IT 



where 

h?A — 2,0-3 ~ 1^16 + |oi4 — Ql8 + 2ai3 
/i^ = a 3 — a 4 + ^(07 — as) + ag + aio + an + ^12 + |(«i6 + «w — a is + a n — 019) 

ft" = — |ai6 + §°14 "I" a 18 "I" 2di 3 

h n A K = a 4 + v^ag + a 9 + \a w - 2a n - a 15 + a u + |(ai 8 + a 19 + a 13 ) (B9) 

Note only aig contributes to i?^ in the parity violating two pions vertices. In the two 
kaons vertices 03-4, 07-3, aio-19 all lead to nonzero contribution to Ra- 



APPENDIX C: A INTERMEDIATE STATES AND EM INSERTIONS 

As noted in Section IV, the amplitudes of Figures 4-6 vanish through 0(1/A£). Below, 
we briefly summarize the reasons behind this result. 



1. PV 7rAiV contribution 



In the case where the A enters as an intermediate state we have the Feynman diagrams 
shown in Figure 4. Since the final and initial states are both nucleons, the two-pion parity 
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violating vertices in Eqs. @2])-([24j) arise first at two- loop order and contribute to the 
nucleon anapole moment at the order of 0(1/A^). Although the PV A^A7r interactions 
nominally contribute at lower order, in this case such contributions vanish up to (9(1/A^). 
The reason is as follows. Each of the parity violating and CP conserving single pion 
vertices has the same Lorentz structure — 275. In the heavy baryon expansion, the relevant 
vertices are obtained by the substitution P + i , y 5 P + , which vanishes. The leading nonzero 
contribution arises at first order in the l/m N expansion. Consequently, its contribution 
to the nucleon anapole moment appears only at 0(l/A^mN), and since in this work we 
truncate at 0(1/A 2 ), the PV tcAN vertices do not contribute. 



2. Magnetic moment insertions 



The nucleon has a large isovector magnetic moment. We thus consider associated 
possible PV chiral loop corrections which lead to a nucleon anapole moment. The relevant 
diagrams are shown in Figure 5. At 0(1/A^) there are only four relevant diagrams Figures 
5a-d. Since the magnetic moment is of 0(1/A X ) and the strong pion baryon vertex is of 
0(1/F n ), the remaining PV vertex must be a Yukawa coupling if the loop is to contribute 
at (9(1/A^) or lower. For the nucleon magnetic moment insertion we have, for example, 

iM 5a + iM 5b = i<f p e li q v v a [Sp, SV 



m^F^ ' J (2tt) d v ■ k v ■ (q + k) k 2 — ml + ie 

(CI) 

where is the photon polarization vector and fi^ is the nucleon magnetic moment. The 
denominator of the integrand in ( |Cl| ) is nearly the same as for M 3e . The numerator 
contains a single factor of S ■ k. Hence, Figures 5a and 5b vanish for the same reason as 
does M 3e . 

For the nucleon delta transition magnetic moment insertion we have 



tM bc + iM 5d = --j= — — (q a e u - e a q u ) [P£ /2 S + S P 3 f 2 ] 

d D k k„ 1 1 



(27r) D v ■ k v ■ (q + k) k 2 — m 2 + ie 



(C2) 



where /iav is the nucleon delta transition magnetic moment and Py 2 = g^— v^v u + ■^S'^S 1 ' 
is the spin | projection operator in the heavy baryon chiral perturbation framework. 
Since the integrand is the same as in M^ e the integral is proportional to t> M . Moreover, 
v^Py 2 = 0, so that M 5c + M 5 d = 0. Finally, the A magnetic insertions of Figures 5e-h 
require the PV A^A7r vertex, which starts off at 0(l/m N F n ). Thus, the latter do not 
contribute up to (9(1/A^). In short, none of the magnetic insertions contribute at the 
order to which we work in this analysis. 
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3. PV electromagnetic insertions 



Another possible source to the nucleon anapole moment arises from the PV magnetic 
moment like insertions as shown in Figure 6. All three PV 'yNN vertices Ci_3 in Eq. 
(pip and PV 7AiV vertices <i 4 _ 6 in Eq. (^) start off with one pion, so they are of order 
0(1/ A x F n ). Vertices d-j^ have two pions and are order 0(1/ k x F%). The corresponding 
contributions to the nucleon anapole moment appear at order of C(l/A^) or C(l/A^), 
respectively. The leading PV '-/AN vertices c?i_3 do not have pions and are of the order 
0(1/ A x ). In our case, however, the final and initial states are both nucleons. The A 
can appear as the intermediate state inside the chiral loop, which leads to an additional 
factor 1/F% from two strong vertices. In the end g?i_3 contributes to the nucleon anapole 
moment at (9(1/A^). Finally, the PV 7AA vertices contain one ir. Since the A can 
only appear as an intermediate state, this vertex contributes at two-loop order and is of 
higher-order in chiral counting than we consider here (the corresponding diagrams are not 
shown). Thus, to 0(1/A X ), the PV electromagnetic insertions do not contribute. 
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Figure Captions 

FIG 1. Axial vector 'jNN coupling, generated by PV hadronic interactions 

Figure 2. Feynman diagrams for polarized electron nucleon scattering. Figure 2a gives 
tree-level Z°-exchange amplitude, while FIG. 2b gives the anapole moment contribution. 
The dark circle indicates an axial vector coupling. 

Figure 3. Meson-nucleon intermediate state contributions to the nucleon anapole 
moment. The shaded circle denotes the PV vertex. The solid, dashed and curly lines 
correspond to the nucleon, pion and photon respectively. For the SU(3) case the 
intermediate states can also be hyperons and kaons. 

Figure 4. The contribution to the nucleon anapole moment from PV 7rAiV vertices. The 

double line is the A intermediate state. 

Figure 5. Anapole moment contributions generated by insertions of the baryon magnetic 
moment operator, denoted by the cross, and the PV hadronic couplings, denoted by the 

shaded circle. 

Figure 6. PV electromagnetic insertions, denoted by the overlapping cross and shaded 

circle. 

Figure 7. Vector meson contribution to the anapole moment. Shaded circle indicates PV 

hadronic coupling. 

Figure 8. Contributions to (a) PV NN interaction and (b) PV two-body current 
generated by the vector terms in Eqs. ( |T8[|20| ). 
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